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Comment on "Origin of Giant Optical Non- 
linearity in Charge- Transfer Mott Insulators: A 
New Paradigm for Nonlinear Optics" 

Zhang jlj has proposed a novel mechanism for the gi- 
ant optical nonlinearities in linear chain Nickel Halides 
(Ni-X) @ and the cuprate Sr 2 Cu0 3 f§. We show that 
Zhang's theory is inapplicable to these systems and also, 
some of his numerical results are finite size artifacts. Al- 
though Zhang's Hamiltonian for Ni-X is correct, the na- 
ture of the eigenstates depend crucially on the relative 
site energies of Ni and X in A ^ .(— Vfn^ a . Zhang- 
chose the Ni atoms to occupy the odd sites in his exact 
diagonalization of N = 8 sites, such that the site energy 
of Ni is lower than that of X. There is then a competition 
between the on-site correlation Um (which prefers Ni 3+ ) 
and A (which prefers Ni 2+ ), and near the interface, the 
nonlinear coefficient 7 is huge for A < A C (U) and zero 
for A > A C (U) (see Fig. 1 of Q). Based on earlier esti- 
mates of U and the magnitude of A, Zhang claims that 
Ni-X lie very close to the interface. 

Zhang's choice of the relative site energies is opposite to 
that of all previous authors and can be precluded. 

Within the Ni 3+ phase, the charge-transfer (CT) gap 
with conventional parametrization QJ^] increases with 
{U Ni — Ux) + V + 2A, whereas with Zhang's choice the 
gap increases with (Um — Ux) + V — 2A. With a less 
diffuse valence atomic orbital in CI, Uci > Ubt- Since, 
however, CI is more electronegative than Br, 2A is larger 
for X = CI with conventional parametrization [QJ^] and 
compensates for the smaller (Um — Ux) to give a larger 
optical gap in Ni-Cl, as is experimentally observed ||. 
With Zhang's choice of site energies, the only way to 
explain the larger CT gap in Ni-Cl is to make the un- 
physical assumption that 2 A is smaller for CI than Br, 
i.e., CI is less electronegative than Br! Thus the con- 
ventional parametrization is correct. Hence there is no 
competition between Um and A, and no enhancement of 

7- 

A second error in |lj] involves the the sudden drop in 7 

by orders of magnitude to zero (see Fig. 1 of jj]). We 
show that this is a consequence of a crossover of the 
ground state from total spin S = to S = 1, a finite 
size effect. It is well known that the ground state of 
a finite undistorted non-half-filled periodic ring with 4n 
electrons is S = 1. With Um 7^ Ux and V > this 
occurs at nonzero A. Zhang must have used a numer- 
ical approach that conserves total S z and not total S. 
In Fig. 1(a) we have plotted AE = E(S Z = 1) - E(S 2 
= 0) for the parameters of Zhang's Fig. 1. In all cases 
AE vanishes at exactly the same A c as in Q. Using a 
method that conserves total S, we have confirmed that 
for A < A c , E(S = 0) = E(S Z = 0), while for A > A c , 
E(S = 0) > E(S Z = 0). With specifically N = 8 and 
12 electrons, the S = 1 ground state is not coupled to 
excited states by the current operator, and this is the 
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FIG. 1. (a) Energy difference between the lowest S z = 1 
and Sz = states of the N = 8 periodic ring with 12 elec- 
trons, for, from left to right, Un% = 2, 3, 4 and 5, with other 
parameters the same as in Fig. 1 of Jl[. (b) The CT gap and 
the matrix element of the current operator J for S = and 
Un% =4. The arrow indicates the A where spin crossover 
occurs. 


reason for the sharp drop in 7 in [jlj . We have confirmed 
that the spin crossover does not occur for N = 12 with 18 
electrons. Also, the true A c that defines the Ni 3+ -Ni 2+ 
interface is larger for all U. In Fig. 1(b) we have plotted 
the CT gap as well as the matrix element of the current 
operator J01 between the lowest S = state and the S = 
one-photon state for the N = 8 periodic ring, using the 
site energies in [jj]. J01 is symmetric about the true A c 
indicating a nearly symmetric behavior of 7 even with 
Zhang's parametrization. 

To conclude: (i) Zhang's choice of site energies for Ni- 
X are incorrect, (ii) his determination of A c is incorrect, 
and (iii) the calculated behavior of 7 for A > A c in [j]] 
is an artifact of finite size calculations. Finally, the giant 
increase in 7 near A C (U) jj| within the model is largely 
due to decreasing CT gap rather than increasing Jqi- 
This does not permit device application, as losses due 
to absorption at wavelengths of interest would be large. 
We acknowledge useful discussions with A. Painelli and 
partial support from NSF-DMR. 
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